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Abstract 

With many Hamiltonians one can naturally associate a l^p-distributed 
Markov process. For nonrelativistic quantum mechanics, this process is in fact 
deterministic, and is known as Bohmian mechanics. For the Hamiltonian of a 
quantum field theory, it is typically a jump process on the configuration space of 
a variable number of particles. We define these processes for regularized quan- 
tum field theories, thereby generalizing previous work of John S. Bell [3] and of 
ourselves We introduce a formula expressing the jump rates in terms of the 
interaction Hamiltonian, and establish a condition for finiteness of the rates. 

PACS numbers: 03.65.Ta (foundations of quantum mechanics), 02.50.-r (proba- 
bility theory, stochastic processes, and statistics), 03.70.+k (theory of quantized 
fields) 
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1 Introduction 

The central formula of this paper is 

, [{2/h)\m{m\P{dq)HP{dq')\n + m 
a[dq \ q) - (*\P(d<f)\*) ' [) 

It plays a role similar to that of Bohm's equation of motion (|2j). Together, these two 
equations make possible a formulation of quantum field theory (QFT) that makes no 
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reference to observers or measurements, while implying that observers, when making 
measurements, will arrive at precisely the results that QFT is known to predict. Special 
cases of formula ((TJ have been utilized before [cO El EH] ■ Part of what we explain in this 
paper is what this formula means, how to arrive at it, when it can be applied, and what 
its consequences are. Such a formulation of QFT takes up ideas from the seminal paper 
of John S. Bell and we will often refer to theories similar to the model suggested by 
Bell in |3, as "Bell-type QFTs". (What similar means here will be fleshed out in the 
course of this paper.) 

The aim of this paper is to define a canonical Bell-type model for more or less any 
regularized QFT. We assume a well-defined Hamiltonian as given; to achieve this, it is 
often necessary to introduce cut-offs. We shall assume this has been done where needed. 
In cases in which one has to choose between several possible position observables, for 
example because of issues related to the Newton- Wigner operator HH] , we shall also 
assume that a choice has been made. 

The primary variables of Bell-type QFTs are the positions of the particles. Bell sug- 
gested a dynamical law, governing the motion of the particles, in which the Hamiltonian 
H and the state vector \1/ determine the jump rates a. We point out how Bell's rates 
fit naturally into a more general scheme summarized by (Q). Since these rates are in a 
sense the smallest choice possible (as explained in Section 03), we call them the minimal 
jump rates. By construction, they preserve the |\l/| 2 distribution. Most of this paper 
concerns the properties and mathematical foundations of minimal jump rates. In Bell- 
type QFTs, which can be regarded as extensions of Bohmian mechanics, the stochastic 
jumps often correspond to the creation and annihilation of particles. We will discuss 
further aspects of Bell- type QFTs and their construction in our forthcoming work [T2j . 

The paper is organized as follows. In Section 2 we introduce all the main ideas 
and reasonings; a superficial reading should focus on this section. Some examples of 
processes defined by minimal jump rates are presented in Section 3. In Section 4 we 
provide conditions for the rigorous existence and finiteness of the minimal jump rates. In 
Section 5 we explain in what sense the rates (0) are minimal. Section 6 concerns further 
properties of processes defined by minimal jump rates. In Section 7 we conclude. 



2 The Jump Rate Formula 

2.1 Review of Bohmian Mechanics and Equivariance 

Bohmian mechanics is a non-relativistic theory about N point particles moving 

in 3-space, according to which the configuration Q = (Q l7 . . . , Q N ) evolves according 



to 1 



-£ = v(Q), v = hlm- W¥ . (2) 



lr The masses rrik of the particles have been absorbed in the Riemann metric g^ v on configuration 
space M. 3N , gi a ,jb = m i $ij $ab, hi — 1...N, a,b = 1, 2, 3, and V is the gradient associated with g M „, 
i.e., V = (m^ 1 V qi ,...,m N 1 V qN ). 
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\I/ = ^t(g) is the wave function, which evolves according to the Schrodinger equation 



t h— = HV, (3) 

with 

H = -^A + V (4) 

for spinless particles, with A = divV. For particles with spin, \1/ takes values in the 
appropriate spin space C fc , V may be matrix valued, and numerator and denominator 
of (j2j) have to be understood as involving inner products in spin space. The secret 
of the success of Bohmian mechanics in yielding the predictions of standard quantum 
mechanics is the fact that the configuration Q t is |\E^| 2 -distributed in configuration space 
at all times t, provided that the initial configuration Q (part of the Cauchy data of the 
theory) is so distributed. This property, called equivariance in suffices for empirical 
agreement between any quantum theory (such as a QFT) and any version thereof with 
additional (often called "hidden") variables Q, provided the outcomes of all experiments 
are registered or recorded in these variables. That is why equivariance will be our guide 
for obtaining the dynamics of the particles. 

The equivariance of Bohmian mechanics follows immediately from comparing the 
continuity equation for a probability distribution p associated with (J2J), 

^ = - div (pv) , (5) 
with the equation satisfied by |\l/| 2 which follows from (|3*|). 

^£ ( g> t ) = \ Im (q, t) (HV) (q, t)] • (6) 
In fact, it follows from (@J) that 



— Im 

n 



V*(q,t)(HV)(q,t) = -div Mm^*(g, t) W(g, t) (7) 
so, recalling (J2J), one obtains that 



dt 



— div (|^r^) , (8) 



and hence that if p t = \^t\ 2 a t some time t then p t = \^t\ 2 for all times. Equivariance 
is an expression of the compatibility between the Schrodinger evolution for the wave 
function and the law, such as (J2J), governing the motion of the actual configuration. In 
[Tlj . in which we were concerned only with the Bohmian dynamics (J2J), we spoke of the 
distribution \^\ 2 as being equivariant. Here we wish to find processes for which we have 
equivariance, and we shall therefore speak of equivariant processes and motions. 
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2.2 Equivariant Markov Processes 

The study of example QFTs like that of fT] has lead us to the consideration of Markov 
processes as candidates for the equivariant motion of the configuration Q for Hamilto- 
nians H more general than those of the form (@J. 

Consider a Markov process Q t on configuration space. The transition probabilities 
are characterized by the backward generator L t , a (time-dependent) linear operator 
acting on functions / on configuration space: 

L t f(q) = ^-E(f(Q t+s )\Q t = q) (9) 

where d/ds means the right derivative at s = and E( ■ | ■ ) denotes the conditional 
expectation. Equivalently, the transition probabilities are characterized by the forward 
generator Jzf t (or, as we shall simply say, generator) , which is also a linear operator but 
acts on (signed) measures on the configuration space. Its defining property is that for 
every process Q t with the given transition probabilities, the distribution p t of Q t evolves 
according to 

< io > 

S£ t is the dual of L t in the sense that 

J f(q) Sftp{dq) = f LJ(q) p(dq) . (11) 

We will use both L t and Sf t , whichever is more convenient. We will encounter several 
examples of generators in the subsequent sections. 

We can easily extend the notion of equivariance from deterministic to Markov pro- 
cesses. Given the Markov transition probabilities, we say that the |\I/| 2 distribution is 
equivariant if and only if for all times t and t' with t < t', a configuration Q t with dis- 
tribution \^t\ 2 evolves, according to the transition probabilities, into a configuration Q t > 
with distribution |^t'| 2 . In this case, we also simply say that the transition probabilities 
are equivariant, without explicitly mentioning l^] 2 . Equivariance is equivalent to 

j?m 2 = ^f (12) 

for all t. When (JT2j) holds (for a fixed t) we also say that ££ t is an equivariant generator 
(with respect to ^ t and H). Note that this definition of equivariance agrees with the 
previous meaning for deterministic processes. 

We call a Markov process Q equivariant if and only if for every t the distribution p t 
of Qt equals \^t\ 2 - For this to be the case, equivariant transition probabilities are nec- 
essary but not sufficient. (While for a Markov process Q to have equivariant transition 
probabilities amounts to the property that if p t = \^t\ 2 for one time t, where p t denotes 
the distribution of Qt-, then p t > = \^t'\ 2 for every t' > t, according to our definition of 
an equivariant Markov process, in fact p t = \^t\ 2 for all t.) However, for equivariant 
transition probabilities there exists a unique equivariant Markov process. 
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The crucial idea for our construction of an equivariant Markov process is to note 
that (jUJ) is completely general, and to find a generator S£ t such that the right hand side 
of (jUJ) can be read as the action of Jzf on p = |\I/| 2 , 

^Imty*H^ = ^\^\ 2 . (13) 

We shall implement this idea beginning in Section |2~H after a review of jump processes 
and some general considerations. But first we shall illustrate the idea with the familiar 
case of Bohmian mechanics. 

For H of the form (0J), we have (|IJ) and hence that 

2 T , T .™. , T ._, ,. A. T ., 2j . T **v^ 



-lm^f*H^f = -div (Mm#*W) = -div M^rMm— - — \ . (14) 

Since the generator of the (deterministic) Markov process corresponding to the dynam- 
ical system dQ/dt = v(Q) given by a velocity vector field v is 

Sfp = - div (pu) , (15) 

we may recognize the last term of (|14|) as ££ |\1/| 2 with «£? the generator of the determin- 
istic process defined by (0). Thus, as is well known, Bohmian mechanics arises as the 
natural equivariant process on configuration space associated with H and \& . 

To be sure, Bohmian mechanics is not the only solution of (fT3|) for H given by 
(J3J). Among the alternatives are Nelson's stochastic mechanics [2H] and other velocity 
formulas jSj. However, Bohmian mechanics is the most natural choice, the one most 
likely to be relevant to physics. (It is, in fact, the canonical choice, in the sense of 
minimal process which we shall explain in ^21 Sec. 5.2].) 

An important class of equivariant Markov processes are equivariant jump processes, 
which we discuss in the next three sections. They arise naturally in QFT, as we shall 
explain in Section 1231 



2.3 Equivariant Jump Processes 

Let Q denote the configuration space of the process, whatever sort of space that may 
be (vector space, lattice, manifold, etc.); mathematically speaking, we need that Q be a 
measurable space. A (pure) jump process is a Markov process on Q for which the only 
motion that occurs is via jumps. Given that Q t = q, the probability for a jump to q', 
i.e., into the infinitesimal volume dq' about q', by time t + dt is cr t (dq'\q) dt, where a is 
called the jump rate. In this notation, o is a finite measure in the first variable; o(B\q) 
is the rate (the probability per unit time) of jumping to somewhere in the set B C Q, 
given that the present location is q. The overall jump rate is cr(Q\q). 

It is often the case that Q is equipped with a distinguished measure, which we shall 
denote by dq or dq', slightly abusing notation. For example, if Q = M. d , dq may be the 
Lebesgue measure, or if Q is a Riemannian manifold, dq may be the Riemannian volume 
element. When a( ■ \q) is absolutely continuous relative to the distinguished measure, 
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we also write cr(q'\q) dq' instead of a(dq'\q). Similarly, we sometimes use the letter p for 
denoting a measure and sometimes the density of a measure, p(dq) = p(q) dq. 

A jump first occurs when a random waiting time T has elapsed, after the time to a t 
which the process was started or at which the most recent previous jump has occurred. 
For purposes of simulating or constructing the process, the destination q' can be chosen 
at the time of jumping, t + T, with probability distribution cr tQ+ r(Q|9) _1 °"to+r( ' !<?)• 
In case the overall jump rate is time-independent, T is exponentially distributed with 
mean o-(Qlq)^ 1 - When the rates are time-dependent — as they will typically be in what 
follows — the waiting time remains such that 



J t 

is exponentially distributed with mean 1, i.e., T becomes exponential after a suitable 
(time-dependent) rescaling of time. For more details about jump processes, see 0. 
The generator of a pure jump process can be expressed in terms of the rates: 



a "balance" or "master" equation expressing dp/dt as the gain due to jumps to dq minus 
the loss due to jumps away from q. 

We shall say that jump rates a are equivariant if ££ a is an equivariant generator. It 
is one of our goals in this paper to describe a general scheme for obtaining equivariant 
jump rates. In Sections 12.41 and 12.51 we will explain how this leads us to formula (0). 

2.4 Integral Operators Correspond to Jump Processes 

What characterizes jump processes versus continuous processes is that some amount 
of probability that vanishes at q 6 Q can reappear in an entirely different region of 
configuration space, say at q' G Q. This is manifest in the equation for dp/dt, (|TB)l: 
the first term in the integrand is the probability increase due to arriving jumps, the 
second the decrease due to departing jumps, and the integration over q' reflects that q' 
can be anywhere in Q. This suggests that Hamiltonians for which the expression © 
for d\^\ 2 /dt is naturally an integral over dq' correspond to pure jump processes. So 
when is the left hand side of ()13|) an integral over dq'l When H is an integral operator, 
i.e., when (q\H\q') is not merely a formal symbol, but represents an integral kernel that 
exists function or a measure and satisfies 



(For the time being, think of Q as M. d and of wave functions as complex valued.) In this 
case, we should choose the jump rates in such a way that, when p = |\1/| 2 , 




^ap{dq) = / (a(dq\q')p(dq') - a(dq'\q)p(dq)) , 



(16) 



q'eQ 




(17) 



a(q\q')p(q')-a(q'\q)p(q) = |lm **(<?) (q\H\<l')W) 



(18) 
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and this suggests, since jump rates must be nonnegative (and the right hand side of (fTHj) 
is anti-symmetric), that 



a(q\q') p(q') 



-2 
Ji 



Im**(?) (q\H\q')m{q) 



+ 



(where x + denotes the positive part of i 6 1, that is, x + is equal to x for x > and is 
zero otherwise), or 



These rates are an instance of what we call the minimal jump rates associated with H 
(and \1/). They are also an instance of formula as will become clear in the following 
section. The name comes from the fact that they are actually the minimal possible 
values given (fTSj) . as is expressed by the inequality (|9l)j) and will be explained in detail 
in Section Minimality entails that at any time t, one of the transitions q\ — > qi or 
<?2 — * <7i is forbidden. We will call the process defined by the minimal jump rates the 
minimal jump process (associated with H). 

In contrast to jump processes, continuous motion, as in Bohmian mechanics, cor- 
responds to such Hamiltonians that the formal matrix elements (q\H\q') are nonzero 
only infinitesimally close to the diagonal, and in particular to differential opera- 
tors like the Schrodinger Hamiltonian (0J), which has matrix elements of the type 



The minimal jump rates as given by ()19|) have some nice features. The possible 
jumps for this process correspond to the nonvanishing matrix elements of H (though, 
depending on the state even some of the jump rates corresponding to nonvanishing 
matrix elements of H might happen to vanish). Moreover, in their dependence on the 
state the jump rates a depend only "locally" upon ^: the jump rate for a given jump 
q' —>■ q depends only on the values ^f(q') and ^(q) corresponding to the configurations 
linked by that jump. Discretizing M 3 to a lattice el? , one can obtain Bohmian mechanics 
as a limit e — > of minimal jump processes whereas greater-than-minimal jump 

rates lead to Nelson's stochastic mechanics j2H| and similar diffusions; see pffi IT7] . If 
the Schrodinger operator (@J) is approximated in other ways by operators corresponding 
to jump processes, e.g., by H E = e~ eH He~ eIi , the minimal jump processes presumably 
also converge to Bohmian mechanics. 

We have reason to believe that there are lots of self-adjoint operators which do not 
correspond to any stochastic process that can be regarded as defined, in any reasonable 
sense, by (119)) . 2 But such operators seem never to occur in QFT. (The Klein-Gordon 
operator \/m 2 c A — h 2 c 2 A does seem to have a process, but it requires a more detailed 
discussion which will be provided in a forthcoming work [T5].) 

2 Consider, for example, H = pcosp where p is the one-dimensional momentum operator —ihd/dq. 
Its formal kernel (q\H\q') is the distribution — f <5'(<? — q' — X) — ^S'(q — q' + X), for which (fTU|) would not 
have a meaning. From a sequence of smooth functions converging to this distribution, one can obtain 
a sequence of jump processes with rates ()19|) : the jumps occur very frequently, and are by amounts of 
approximately ±1. A limiting process, however, does not exist. 



[{2/h) Im**(g) (q\H\q')^(q')] 



+ 




(19) 



6' 



(q-q')+V(q)6(q-q f ). 
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2.5 Minimal Jump Rates 



The reasoning of the previous section applies to a far more general setting than just con- 
sidered: to arbitrary configuration spaces Q and "generalized observables" — POVMs — 
defining, for our purposes, what the "position representation" is. We now present this 
more general reasoning, which leads to formula (fT|). 

The process we construct relies on the following ingredients from QFT: 

1. A Hilbert space J$? with scalar product (^|$). 

2. A unitary one-parameter group U in J%f with Hamiltonian H, 

so that in the Schrodinger picture the state \& evolves according to 

*% = H* t . (20) 

U t could be part of a representation of the Poincare group. 

3. A positive-operator-valued measure (POVM) P(dq) on Q acting on Jjf, so that 
the probability that the system in the state is localized in dq at time t is 

F t {dq) = (%\P{dq)\%) . (21) 



Mathematically, a POVM P on Q is a countably additive set function ("measure"), 
defined on measurable subsets of Q, with values in the positive (bounded self-adjoint) 
operators on (a Hilbert space) Jrff, such that P{Q) is the identity operator. 3 Physically, 
for our purposes, P(-) represents the (generalized) position observable, with values in 
Q. The notion of POVM generalizes the more familiar situation of observables given 
by a set of commuting self-adjoint operators, corresponding, by means of the spectral 
theorem, to a projection- valued measure (PVM): the case where the positive operators 
are projection operators. A typical example is the single Dirac particle: the position 
operators on L 2 (M 3 ,C 4 ) induce there a natural PVM Po(-)'- for any Borel set B C 
1R 3 , Pq(B>) is the projection to the subspace of functions that vanish outside B, or, 
equivalently, P (B)^/(q) = l B (q)^(q) with 1 B the indicator function of the set B. 
Thus, (ty\PQ(dq)\\f/) = \^(q)\ 2 dq. When one considers as Hilbert space Jif only the 
subspace of positive energy states, however, the localization probability is given by 
P( ■ ) = p + p ( . )j with P+ : L 2 (R 3 , C 4 ) -> the projection and / : -> L 2 (R 3 , C 4 ) 
the inclusion mapping. Since P + does not commute with most of the operators Pq(B), 
P( ■ ) is no longer a PVM but a genuine POVM 4 and consequently does not correspond to 

3 The countable additivity is to be understood as in the sense of the weak operator topology. This 
in fact implies that countable additivity also holds in the strong topology. 

4 This situation is indeed more general than it may seem. By a theorem of Naimark p. 142], every 
POVM P{ ■ ) acting on J? is of the form P( ■ ) = P+P { ■ )P+ where P is a PVM on a larger Hilbert 
space, and P+ the projection to Jf°. 
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any position operator — although it remains true (for \l/ in the positive energy subspace) 
that ($f\P(dq)\ty) = \^(q)\ 2 dq. That is why in QFT, the position observable is indeed 
more often a POVM than a PVM. POVMs are also relevant to photons |T[ |22]- in one 
approach, the photon wave function \1/ : R 3 — > C 3 is subject to the constraint condition 
V • * = <9i^i + d 2 ^2 + d 3 ^ 3 = 0. Thus, the physical Hilbert space J? is the (closure 
of the) subspace of L 2 (IR 3 ,C 3 ) defined by this constraint, and the natural PVM on 
L 2 (R 3 ,C 3 ) gives rise, by projection, to a POVM on Jflf. So much for POVMs. Let us 
get back to the construction of a jump process. 

The goal is to specify equivariant jump rates o = a^ ,H,p ', i.e., such rates that 

dP 

^P=-. (22) 

To this end, one may take the following steps: 
1. Note that 

dFt{dq) -hm(%\P(dq)H\%). (23) 



dt h 

2. Insert the resolution of the identity / = f P(dq') and obtain 

q'eQ 



dF t (dq) 



dt 

q'eQ 



t (dq,dq'), (24) 



where 

Udq,dq') = \ Im {%\P{dq)HP{dq')\%) . (25) 
n 

3. Observe that J is anti-symmetric, §(dq',dq) = —§(dq,dq f ). Thus, since x = x + — 

(-*) + , 

I(dq,dq') = [(2/ft) Im (^\P(dq)HP(dq')\^}} + - [(2/H) lm(^\P(dq')HP(dq)\^}} + . 

4. Multiply and divide both terms by P( • ), obtaining that 

q'eQ q'eQ 

l(2/h)Im(y\P(dq')HP(dq)m + m ^^ 
(*|P(dg)|*> [q) 

5. By comparison with recognize the right hand side of the above equation as 
Jzf CT P, with S£ a the generator of a Markov jump process with jump rates which 
we call the minimal jump rates. We repeat the formula for convenience: 

a{dq \ q i) = [m)^^\P(dq)HP(dq'm] + 



(V\P(dq')\V) 
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Mathematically, the right hand side of this formula as a function of q' must be understood 
as a density (Radon-Nikodym derivative) of one measure relative to another. The plus 
symbol denotes the positive part of a signed measure; it can also be understood as 
applying the plus function, x + = max(x, 0), to the density, if it exists, of the numerator. 

To sum up, we have argued that with H and \I/ is naturally associated a Markov 
jump process Q t whose marginal distributions coincide at all times by construction with 
the quantum probability measure, p t ( ■ ) = ¥ t ( • ), so that Q t is an equivariant Markov 
process. 

In Section HJ we establish precise conditions on H, P, and \I/ under which the jump 
rates (PJ) are well-defined and finite P- almost everywhere, and prove that in this case 
the rates are equivariant, as suggested by the steps 1-5 above. It is perhaps worth 
remarking at this point that any H can be approximated by Hamiltonians H n (namely 
Hilbert-Schmidt operators) for which the rates (Q) are always (for all \l/) well-defined 
and equivariant, as we shall prove in Section f4. 2. II 



2.6 Bell-Type QFT 

A Bell-type QFT is about particles moving in physical 3-space; their number and posi- 
tions are represented by a point Q t in configuration space Q, with Q defined as follows. 
Let TM. 3 denote the configuration space of a variable (but finite) number of identical 
particles in R 3 , i.e., the union of (R 3 )™ modulo permutations, 

oo 

m 3 = |J(M 3 )7S„. (26) 

n=0 

Q is the Cartesian product of several copies of TIR 3 , one for each species of particles. 
For a discussion of the space TR 3 , and indeed of TS for any other measurable space 5* 
playing the role of physical space, see O Sec. 2.8]. 

A related space, for which we write T^R 3 , is the space of all finite subsets of R 3 ; 
it is contained in TR 3 , after obvious identifications. In fact, T^R 3 = TR 3 \ A, where 
A is the set of coincidence configurations, i.e., those having two or more particles at 
the same position. T^R 3 is the union of the spaces for n = 0, 1, 2, . . ., where 
is the space of subsets of R 3 with n elements, a manifold of dimension 3n (see ^U] 
for a discussion of Bohmian mechanics on this manifold). The set A of coincidence 
configurations has codimension 3 and thus can usually be ignored. We can thus replace 
TR 3 by the somewhat simpler space T^R 3 . 

Qt follows a Markov process in Q, which is governed by a state vector \l/ in a suitable 
Hilbert space Jif. is related to Q by means of a PVM or POVM P. 

The Hamiltonian of a QFT usually comes as a sum, such as 

H = H + Hj (27) 

with H the free Hamiltonian and Hj the interaction Hamiltonian. If several particle 
species are involved, H is itself a sum containing one free Hamiltonian for each species. 
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The left hand side of ()13j) . which should govern our choice of the generator, is then also 
a sum, 

-Imty*H ^ + -Im^H^ = ££W . (28) 

h h 

This opens the possibility of finding a generator 5£ by setting Jzf = Jzf + provided 
we have generators Jzfo arid S£i corresponding to Hq and Hi in the sense that 

-2b M 2 (29a) 
Jzf 7 |^| 2 . (29b) 

This feature of (fT^j) we call process additivity; it is based on the fact that the left hand 
side of f!13|) is linear in H. 

In a Bell-type QFT, the generator J?f is of the form Jzf = J2f +=2/, where Jzfo is usually 
the generator of a deterministic process, usually defined by the Bohmian or Bohm-Dirac 
law of motion, see below, and Jzfj is the generator of a pure jump process, which is our 
main focus in this paper. The process generated by ££ is then given by deterministic 
motion determined by Jz?o, randomly interrupted by jumps at a rate determined by J2?/. 

We thus need to define two equivariant processes, one (the "free process" ) associated 
with Hq and the other (the "interaction process") with Hj. The interaction process is 
the pure jump process with rates given by (HJ) with Hi in place of H. We now give a 
description of the free process for the two most relevant free Hamiltonians: the second- 
quantized Schrodinger operator and the second-quantized Dirac operator. We give a 
more general and more detailed discussion of free processes in [T2]; there we provide a 
formula, roughly analogous to (jl]). for J2f in terms of Hq, and an algorithm for obtaining 
the free process from a one-particle process that is roughly analogous to the "second 
quantization" procedure for obtaining Hq from a one-particle Hamiltonian. 

The free process associated with a second-quantized Schrodinger operator arises from 
Bohmian mechanics. Fock space Jrf? = & is a direct sum 

oo 

jr = 0jr("), (30) 

n=0 

where is the n-particle Hilbert space. J^"-) is the subspace of symmetric (for 

bosons) or anti-symmetric (for fermions) functions in L 2 (R 3n , (C 2s + 1 )« m ) for spin-s par- 
ticles. Thus, \& G & can be decomposed into a sequence ^ = (^(°\ v&W, . . . , vl/w, . . .), 
the n-th member being an n-particle wave function, the wave function represent- 
ing the n-particle sector of the quantum state vector. The obvious way to obtain a 
process on Q = TM. 3 is to let the configuration Q(t), containing iV = #Q(t) particles, 
move according to the iV-particle version of Bohm's law (J2J), guided by ty( N \ 5 This is 

5 As defined, configurations are unordered, whereas we have written Bohm's law for ordered 
configurations. Thanks to the (anti-)symmetry of the wave function, however, all orderings will lead to 
the same particle motion. For more about such considerations, see our forthcoming work 



-lmty*H ty = 
n 

-Im^H^ = 

h 
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indeed an equivariant process since H has a block diagonal form with respect to the 
decomposition (J3DJ), 

oo 
n=0 

and Hq^ is just a Schrodinger operator for n noninteracting particles, for which, as we 
already know, Bohmian mechanics is equivariant. We used a very similar process in ^JJ 
(the only difference being that particles were numbered in [TT]). 

Similarly, if Ho is the second quantized Dirac operator, we let a configuration Q with 
N particles move according to the usual iV-particle Bohm-Dirac law 01 p. 274] 

dQ = vT(Q) a N *(Q) 

dt ° tt*(Q)tt(Q) 1 ' 

where c denotes the speed of light and = (c^ 1 ), . . . , ct^) with a^ k ' acting on the 
spin index of the k-th particle. 

This completes the construction of the Bell-type QFT. An explicit example of a 
Bell-type process for a simple QFT is described in [IT], which we take up again in 
Section f3.12l below to point out how its jump rates fit into the scheme (JJJ). Another 
such example, concerning electron-positron pair creation in an external electromagnetic 
field, is described in [J21 Sec. 3.3.]. 



3 Examples 

In this section, we present various special cases of the jump rate formula (JJJ) and examples 
of its application. We also point out how the jump rates of the models in 11] and J3] 
are contained in (JJJ). 



3.1 A First Example 

To begin with, we consider Q = R d , = L 2 (R d , C), and P the natural PVM, which 
may be written P(dq) = \q)(q\dq. Then, F(dq) = (^\P(dq)\^) = \^(q)\ 2 dq, and the 
jump rate formula (JJJ) reads 



a(q\q') 



[(2/h) Im#*(g) (q\H\q')^(q')]- 
•q/*(q>) ^j(q') 

2 Im tt*(g) (gjgjgQ 
h ty*(q') 



(32a) 
(32b) 



Note that ()32a|) is the same expression as ()19jl. As a simple example of an operator Hj 
with a kernel, consider a convolution operator, Hi = V*, where V may be complex- 
valued and V(— q) = V*(q), 

(Hr*)(q)= [ V{q-q')^{q')dq'. 
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The kernel is (q\Hj\q') = V(q — q'). Together with H = — 4-A, we obtain a baby 
example of a Hamiltonian H = H + Hj that goes beyond the form (@J) of Schrodinger 
operators, in particular in that it is no longer local in configuration space. Recall that 
Hq is associated with the Bohmian motion (J2J). Combining the two generators on the 
basis of process additivity, we obtain a process that is piecewise deterministic, with 
jump rates (|T^|) and Bohmian trajectories between successive jumps. 



3.2 Wave Functions with Spin 

Let us next become a bit more general and consider wave functions with spin, i.e., 
with values in C k . We have Q = R d ,Jf = L 2 (R d , C k ) and P the natural PVM, which 
may be written P{dq) = Y^h=i \q^){Q^\ dq, where i indexes the standard basis of C fc . 
Another way of viewing P is to understand Jf? as the tensor product L 2 (M. d , C) <8> C fc , 
and P(dq) = P (dq) ® I & with P the natural PVM on L 2 (R d , C) and I c * the identity 
operator on C h . Using the notation ((<&(q)\ty(q))) for the scalar product in C fc , we can 
write F(dq) = (^\P(dq)\^/) = (q)\ty (q))) dq, and the jump rate formula (JTJ) reads 

, nlfA _ [(2/h)lm(^(q)\K(q,q')\^(q')W 

{qlq ] ~ «WS (33) 

with K(q,q'), the kernel of H, a k x matrix. If we write $*((?) *&{q) for 
as we did in (J2J) and (|3T|). and for K(q,q'), f)33j) reads 

[(2/fi)Imtt*( g ) ( g |#|</> 



which is ()19|) again, interpreted in a different way. 



3.3 Vector Bundles 

Next consider, instead of the fixed value space C k , a vector bundle E over a Riemannian 
manifold Q, and cross-sections of E as wave functions. In order to have a scalar product 
of wave functions, we need that every bundle fiber E q be equipped with a Hermitian 
inner product (( • | • )),. We consider = L 2 (E) (the space of square-integrable cross- 
sections) and P the natural PVM. For any q and q', K(q, q') then has to be a C-linear 
mapping E q i — > E q , so that the kernel of if is a cross-section of the bundle |J , E q ®E*, 
over Q x Q. (JTJ) then reads 

, I /\ = [(2/fe)Im(^( g )|^( g ^)^(^))) g ] + 

In the following we will use the notation <&*(q) ^{q) for (($>(q)\ty(q))) q and for 
K(q, q'), so that 

[(2/ft)Imtt*(g) ( g |g|g^( g Q] + 

which looks like (J19)) again. 
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3.4 Kernels of the Measure Type 

The kernel (q\H\q') can be less regular than a function. Since the numerator of is 
a measure in q and q', the formula still makes sense (for P the natural PVM) when 
the kernel (q\H\q') is a complex measure in q and q'. The mathematical details will 
be discussed in Section 14.21 For instance, the kernel can have singularities like a Dirac 
5, but it cannot have singularities worse than 5, such as derivatives of 5 (as would 
arise from an operator whose position representation is a differential operator). It can 
happen that the kernel is not a function but a measure even for a very well-behaved 
(even bounded) operator. For example, this is the case for H a multiplication operator 
(i.e., a function V(q) of the position operator), (q\H\q r ) = V(q) 5(q — q'). Note, though, 
that multiplication operators correspond to zero jump rates. 

A nontrivial example of an operator with 5 singularities in the kernel is H = 1 — 
cos(p/po) where p = —ihd/dq is the momentum operator in one dimension, M 3 = 
L 2 (R, C), and p is a constant. The dispersion relation E = 1 — cos(p/p ) begins at p = 
like \{j>/pq) 2 but deviates from the parabola for large p. In the position representation, 
H is the convolution with ((27r) _1/ ' 2 times) the inverse Fourier transform of the function 
1 — cos(hk/p ), and thus (q\H\q') = S(q — q') — ~ 5{q — q' + -^) — ~ 8{q — q' — ^). In this 
case, (JTJ) leads to 

•<*)- Jfc ^^ t ('(.-^i) + '(.^-i))- w 

(Note that nonnegative factors can be drawn out of the plus function.) This formula 
may be viewed as contained in (JT§|| as well, in a formal sense. As a consequence of (J23J), 
only jumps by an amount of ±^ can occur in this case. 

3.5 Infinite Rates 

There also exist Markov processes that perform infinitely many jumps in every finite 
time interval (e.g., Glauber dynamics for infinitely many spins). These processes, which 
we do not count among the jump processes, may appear pathological, and we will not 
investigate them in this paper, but we note that some Hamiltonians may correspond to 
such processes. They could arise from jump rates a{ ■ \q') given by (JIJ) that form not a 
finite but merely a cr-finite measure, so that cr(Q\q') = oo. Here is an (artificial) example 
of cr-finite (but not finite) rates, arising from an operator H that is even bounded. 

Let Q = R, = L 2 (R) with P( ■ ) the position PVM, and let H, in Fourier 
representation, be multiplication by f(k) = a/vt/2 sign(fc). H is bounded since / is. / 
is the Fourier transform of i/x, understood as the distribution defined by the principal 
value integral. As a consequence, H has, in position representation, the kernel (q\H\q f ) = 
i/(q — q'). From (JT§j) we obtain the jump rates 

< I '\ 2 1 l Re^*(q)^(q') ] + 
aiqlq)= h **(q>)*(q>) l q-q> \ ' (36) 
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which entails that cr(R|g') = J <r(q\q') dq = oo at least whenever \1/ is continuous (and 
nonvanishing) at q'. Nonetheless, since the rate for jumping anywhere outside the inter- 
val \q' — e,q' + e] is finite for every e > and since J" g f_ £ s \q — q'W{q\q') dq < oo, a process 
with these rates should exist: among the jumps that the process would have to make 
per unit time, the large ones would be few and the frequent ones would be tiny — too 
tiny to significantly contribute. 



3.6 Discrete Configuration Space 

Now consider a discrete configuration space Q. Mathematically, this means Q is a 
countable set. In this case, measures are determined by their values on singletons {q}, 
and we can specify all jump rates by specifying the rate cr(q\q') for each transition q' — > q. 
(Q) then reads 

[(2/h) Im (y\P{q}HP{q'}\^) 



a(q\q') 



(37) 



<*|PMI*> 

We begin with the particularly simple case that there is an orthonormal basis of 
labeled by Q, {\q) : q G Q}, and P is the PVM corresponding to this basis, P{q} = 
\q)(q\- In this case, the notation (q\H\q') and the name "matrix element" can be taken 
literally. The rates (0) then simplify to 



a(q[q') 



[(2/h) Im($\q){q\H\q')(q'\y))- 

WWW) 
2 (*\q)(q\H\</) 



(38a) 
(38b) 



Note that (|38aJ) is the obvious discrete analogue of (fT9|) : in fact, one can regard (|T9|) as 
another way of writing ()38aj) in this case. 

Consider now the more general case that a basis of Hilbert space is indexed by two 
"quantum numbers," the configuration q and another index i. Then the POVM is given 
by the PVM P{q} = J2i l<2S the projection onto the subspace associated with q 

(whose dimension might depend on q); such a PVM may be called "degenerate." We 
have F(q) = (^\P{q}\^) = ^(*|g,i)(g,i|*), and (JTJ) becomes 



a(q\q') 



lmJ2(*\q,i)(q,i\H\q',i f )(q f ,i f \V) 



(39) 



We may also write ((HS)) as (|38ajl . understanding (&\q) and (q'\^) as multi-component, 
(q\H\q') as a matrix, and products as inner products. In case that the dimension of the 
subspace associated with q is always k, independent of q, (|3~9"j) is a discrete analogue of 
the rate formula (}33|) for spinor-valued wave functions. 

Apart from serving as mathematical examples, discrete configuration spaces are rel- 
evant for several reasons: First, they provide particularly simple cases of jump processes 
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with minimal rates that are easy to study. Second, any numerical computation is dis- 
crete by nature. Third, one may consider approximating or replacing the M 3 that is 
supposed to model physical space by a lattice Z 3 ; after all, lattice approaches have often 
been employed in QFT, for various reasons. Moreover, Bell-type QFTs will usually have 
as configurations the positions of a variable number of particles; so the configuration 
has a certain continuous aspect, the positions, and a certain discrete aspect, the number 
of particles. Sometimes one wishes to study simplified models, and in this vein it may 
be interesting to have only the particle number as a state variable, and thus the set of 
nonnegative integers as configuration space. 

3.7 Bell's Process 

The model Bell specified in ^ is a case of a minimal jump process on a discrete set. "For 
simplicity," Bell considers a lattice A instead of continuous 3-space, and a Hamiltonian 
of a lattice QFT. As a consequence, the configuration space Q = T(A) is countable. 
(Bell even makes Q finite, but this is not relevant here. We also remark that according 
to Bell's formulation, even distinguishable particles have configuration space r(A).) 

Bell chooses as the configuration the number of fermions at every lattice site, rather 
than the total particle number (i.e., in our terminology he takes P{q} to be the pro- 
jection to the joint eigenspace of the fermion number operators for all lattice sites with 
eigenvalues the occupation numbers corresponding to q G T(A)). He thus gives the 
fermionic degrees of freedom a status different from the bosonic ones. That is to say, 
boson particles do not exist in Bell's model, despite the fact that Jif = fermions® ^bosons 
and the presence of bosonic terms in the Hamiltonian. 

Thus the PVM P{q} = Pf C rmions{<?} <%> lbosons is "doubly" degenerate: the fermionic 
occupation number operators do not form a complete set of commuting operators, be- 
cause of both the spin and the bosonic degrees of freedom. Different spin states and 
different quantum states of the bosonic fields are compatible with the same fermion 
occupation numbers. So a further index i is necessary to label a basis {|g, i)} of J^. 
The jump rates Bell prescribes are then (jH^j) . and are thus a special case of (0). We 
emphasize that here the index i does not merely label different spin states, but states 
of the quantized radiation as well. 

3.8 A Case of POVM 

Consider for Jtf the space of Dirac wave functions of positive energy. The POVM P( • ) 
we defined on it in Section 1231 is. as we have already remarked, not a PVM but a genuine 
POVM and arises from the natural PVM P ( • ) on L 2 (M 3 , C 4 ) by P( • ) = P + P ( • )/ with 
P + : L 2 (R 3 ,C 4 ) -> J$? the projection and / : Jtf -> L 2 (IR 3 ,C 4 ) the inclusion mapping. 
We can extend any given interaction Hamiltonian H on ^ to an operator on L 2 (K 3 , C 4 ), 
P C xt = IHP+. If P cx t possesses a kernel (q\H ext \q'), then H corresponds to a jump 
process, and the rates (JI} can be expressed in terms of this kernel, since for ^ G Jff, 
(*\P(dq)HP(dq')\V) = (y\P + P (dq)IHP + P (dq')I\y) = {^\P (dq)H ext P (dq')\^) = 
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\&*(g) (q\H ext \q') \I>(g') dqdq'. We thus obtain 

, , ,x [(2/ft)Img*(g) (gRx^O^gO, 

a(g|g) = ^pa • (40) 

This POVM is used in the pair creation model of \Y2\ Sec. 3.3]. 



3.9 Another Case of POVM 

Let 3^ = L 2 (R d ) and let P ( • ) be the natural PVM. We obtain a POVM P by smearing 
out Pq with a profile function ip : M d — > [0, oo) with f p(q) dq = 1 and <p(—q) = <p(q), 
e.g., a Gaussian: 

P(5)= y J <p(q'-q)P (dq'). (41) 



9G-B g ' 

Whereas Pq(B) is multiplication by -P(-B) is multiplication by p * 1 B . It leads to 
F(dq) = (<p-k\V\ 2 )(q) dq. 

The jump rate formula (JTJ then yields 

, , /x [(2/fi) Im ( dq" f dq'" cp(q" - q)V*(q") (q"\H\q"') ty(q"') cp(q"' - q')} + 
(710 ) = — — 

vy|y ; / dg" y?(g" - q>) **(g") #(g") 

i.e., the denominator gets smeared out with ip, and the square bracket in the numerator 
gets smeared out with p in each variable. 



3.10 Identical Particles 

The n-particle sector of the configuration space (without coincidence configurations) 
of identical particles T^(R 3 ) is the manifold of n-point subsets of M 3 ; let Q be this 
manifold. The most common way of describing the quantum state of n fermions is by 
an anti-symmetric (square-integrable) wave function f on Q := M 3n ; let Jff be the 
space of such functions. Whereas for bosons $ could be viewed as a function on Q, for 
fermions \& is not a function on Q. 

Nonetheless, the configuration observable still corresponds to a PVM P on Q: for 
B C Q, we set P(B)V(q 1} ...,q n ) = V(q x , ...,q n ) if {<?i, ...,<?„} G 5 and zero oth- 
erwise. In other words, P(B) is multiplication by the indicator function of 7r~ 1 (i?) 
where ir is the obvious projection mapping Q \ A — > Q, with A the set of coincidence 
configurations. 

To obtain other useful expressions for this PVM, we introduce the formal kets |g) 
for q G Q (to be treated like elements of L 2 (Q)), the anti-symmetrization operator S 
(i.e., the projection L 2 (Q) — > Jt?), the normalized anti-symmetrizer 6 s = VnJ.S, and 

6 The name means this: since S is a projection, Svl/ is usually not a unit vector when ^ is. Whenever 
"J G L 2 (Q) is supported by a fundamental domain of the permutation group, i.e., by a set fl C Q on 
which (the restriction of) 7r is a bijection to Q, the norm of S\l/ is l/^fw\, so that is again a unit 
vector. 
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the formal kets \sq) := s\q) (to be treated like elements of J^). The \q) and \sq) are 
normalized in the sense that 

(qW) = S(q ~ <H and (sq\sq') = (-1)^') 5(q - q'), 

where q = vr(g), q' = vr(g'), g(q,q') is the permutation that carries q into q' given that 
q = q', and (— l) e is the sign of the permutation g. Now we can write 



P(dq) 



E 

<?e7r !(g) 



o?g = n! c?g = \sq){sq\ dq, 



(42) 



where the sum is over the n! ways of numbering the n points in q; the last two terms 
actually do not depend on the choice of q G 7T _1 (g), the numbering of q. 
The probability distribution arising from this PVM is 



\^(q)\ 2 dq = n\ \^(q)\ 2 dq = \(sq\^)\ 2 dq 



with arbitrary q G ti 1 (q). 

If an operator H on L 2 (Q) is permutation invariant, 



(43) 



U 1 HU g = H for every permutation g, 



(44) 



where U s is the unitary operator on L 2 (Q) performing the permutation g, then H maps 
anti-symmetric functions to anti-symmetric functions, and thus defines an operator H 
on M'. If H has a kernel (q\H\q') then the kernel is permutation invariant in the sense 
that 

(g(q)\H\g(q')) = (q\H\q') Vf?, (45) 
where g(q X) . . . , q n ) := (<7 e (i), • • • , Q e ( n ))i an d H a l so possesses a kernel, 

(sq\H\sq') = n\ (q\SHS\q') = 1 V(£(g)|tf |f/(<f )>• 

n! z — ' 

In this case O yields 



a(q\q') 



(#|s<f)(s(f I*) 



(46a) 



(46b) 



where q G 7r x (g) and <f G 7r ^(j'), as running variables in ()46a|) and as arbitrary but 
fixed in (l46bl . 
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3.11 Another View of Fermions 

There is a way of viewing fermion wave functions as being defined on Q, rather than 
M 3n , by regarding them as cross-sections of a particular 1-dimensional vector bundle 
over Q. To this end, define an n!-dimensional vector bundle E by 

E q := C. (47) 

Every function \1/ : M 3n — > C naturally gives rise to a cross-section $ of E, defined by 

*(?):= *(?)■ (48) 

The anti-symmetric functions form a 1-dimensional subbundle of E (see also ^U] for 
a discussion of this bundle). The jump rate formula for vector bundles fJMj) can be 
applied to either the subbundle or E, depending on the way in which the kernel of H is 
given. The kernel (q\H\q') above translates directly into a kernel onQxQ with values 
in E q <S) E*,, for which the rate formula for bundles (|34j) is the same as the rate formula 
for identical particles (J46a|) derived in the previous section. 

Another alternative view of a fermion wave function is to regard it as a complex 
differential form of full rank, a 3n-form, on Q. (See, e.g., jlOj . This would not work 
if the dimension of physical space were even.) Of course, the complex 3n-forms are 
nothing but the sections of a certain 1-dimensional bundle, usually denoted C <8> A 3n <2, 
which is equivalent to the subbundle of E considered in the previous paragraph, and 
which is contained in the bundle C ® AQ of Grassmann numbers over Q. 

3.12 A Simple QFT 

We presented a simple example of a Bell-type QFT in ^Tj , and we will now briefly point 
to the aspects of this model that are relevant here. The model is based on one of the 
simplest possible QFTs jUl p. 339]. 

The relevant configuration space Q for a QFT (with a single particle species) is 
the configuration space of a variable number of identical particles in M 3 , which is the 
set T(]R 3 ), or, ignoring the coincidence configurations (as they are exceptions), the set 
r^(IR 3 ) of all finite subsets of R 3 . The n-particle sector of this is a manifold of dimension 
3n; this configuration space is thus a union of (disjoint) manifolds of different dimensions. 
The relevant configuration space for a theory with several particle species is the Cartesian 
product of several copies of r^(lR 3 ). In the model of [TT], there are two particle species, 
a fermion and a boson, and thus the configuration space is 

Q = T^M 3 ) x r^(M 3 ). (49) 

We will denote configurations by q = (x, y) with x the configuration of the fermions and 
y the configuration of the bosons. 
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For simplicity, we replaced in JT] the sectors of r^(IR 3 ) xT^(IR 3 ), which are manifolds, 
by vector spaces of the same dimension (by artificially numbering the particles), and 
obtained the union 

oo oo 

Q= \J(R 3 ) n x |J(M 3 )"\ (50) 

n=0 m=0 

with n the number of fermions and m the number of bosons. Here, however, we will 
use as the configuration space. In comparison with (J5U|) . this amounts to (merely) 
ignoring the numbering of the particles. 

ffl is the tensor product of a fermion Fock space and a boson Fock space, and thus the 
subspace of wave functions in L 2 (Q) that are anti-symmetric in the fermion coordinates 
and symmetric in the boson coordinates. Let S denote the appropriate symmetrization 
operator, i.e., the projection operator L 2 (Q) — > and s the normalized symmetrizer 

s*(xi, ...,x n ,y 1 ,...,y m ) = Vn\m\S^(x 1} . . . , x n , y x , . . . , y m ), (51) 

i.e., s = V N\ M\ S with N and M the fermion and boson number operators, which 
commute with S and with each other. As in Section TJ. 101 we denote by 7r the projection 
mapping Q \ A -> Q, n(x u . . . , x n , y x , . . . , y m ) = ({x u . . . , x n }, {y x , . . . , y m }). The 
configuration PVM P(B) on Q is multiplication by l^-i^), which can be understood as 
acting on though it is defined on L 2 (Q), since it is permutation invariant and thus 
maps J$? to itself. We utilize again the formal kets \q) where q G Q \ A is a numbered 
configuration, for which we also write q = (x, y) = (x%, . . . , x n , y X) . . . , y m ). We also use 
the symmetrized and normalized kets \sq) = s\q). As in (j42)l . we can write 

P(dq) = \q)(q\dq = n\m\ S\q)(q\dq = \sq)(sq\dq (52) 

with arbitrary q G 7r _1 (g). For the probability distribution, we thus have, as in (J43|) . 

F(dq) = \y(q)\ 2 dq = n\m\\ty(q)\ 2 dq=\(sq\^)\ 2 dq (53) 

with arbitrary q G 7r _1 (g). 

The free Hamiltonian is the second quantized Schrodinger operator (with zero poten- 
tial), associated with the free process described in Section l2~o1 The interaction Hamil- 
tonian is defined by 

Hj = J d 3 xij\x)(al(x) + a v (x))ij(x) (54) 

with ip\x) the creation operators (in position representation), acting on the fermion 
Fock space, and a^(x) the creation operators (in position representation), acting on the 
boson Fock space, regularized through convolution with an L 2 function (p : M. 3 — > R. Hi 
has a kernel; we will now obtain a formula for it, see (jOOj) below. The \sq) are connected 
to the creation operators according to 

\sq) = tf{x n ) ■ ■ -^(xja^yj ■ ■■a\y 1 )\0) , (55) 
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where |0) G J^ 5 denotes the vacuum state. A relevant fact is that the creation and annihi- 
lation operators ip\ij),(tf and a possess kernels. Using the canonical ( ant i-) commutation 
relations for ip and a, one obtains from (|55|) the following formulas for the kernels of 
ip(r) and a(r), r G M 3 : 

(sq\iP(r)\sq') = 5 n>n ^ 8 m , m , 5 3n ' (x Ur-x') 5 ^ y _ yl) (56) 

(sq\a(r)\sq') = 6 n , n , 6 3n (x - x') (-1)^') 6 3m ' (y Ur-y') (57) 

where (x, y) = q = vr(g), and g(x, x') denotes the permutation that carries x to x' given 
that x = x' . The corresponding formulas for ip^ and can be obtained by exchanging 
q and q' on the right hand sides of (JoT)|) and (joT|) . For the smeared-out operator a^(r), 
we obtain 

(sqK(r)\sq>) = 6 n , n , 6 m , m >-i 5 3n (x - x') (-1)^') ^ 5 3 ™( y - y' \ y>) tp(y> - r) (58) 
We make use of the resolution of the identity 



dq \sq)(sq\ 



(59) 



s 



Inserting (p)9"j) twice into (jB^j) and exploiting (|5T)|) and (|55)l. we find 

(sqlHjlsq'} = 5 n , n > <5 m _ w 6 3n (x - x') (-1)^*') ^ 5 3m '(y \ y - y') E ^ - x ) 



+ $n,n> Sm>-l,m ^(x - *') (-1)^ ]T 5 3m (y - y' \ y') £ ^(y' - x) . 

(60) 



This is another case of a kernel containing 5 functions (see Section 
By (J52J), the jump rates (JTJ) are 



a{q\q>) 



|lm (sg|i7j|sg') (s^'l^) 



{^\sq')(sq'\ty) 
More explicitly, we obtain from (J6Uj) the rates 

<r(qW) = <W $ m -i, m > 5 3n (x - x')^5 3m> {y \y-y') a cxea ,{q' U y\q') 



(61) 



+ <W W'-i 5 3n (x - x') £ 5 3m (y - y' \ y') a ann (q' \ y'\q') 



(62) 



with 



Crcrea(g' U y|# 
^ann(g / \ y'\q 



Im**(g) (-l)f(^') £ <^(y - a;') *(g') 



/J. 



^*(g') 



Im^*(g) (-1)^') Y, V(y' -x')V(q') 



(63a) 
(63b) 
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for arbitrary q' G 7r _1 (g') and q G vr _1 (g) with q = (x', y'Uy) respectively q = (x', y'\y'). 
(Note that a sum sign can be drawn out of the plus function if the terms have disjoint 
supports.) 

Equation (|62|) is worth looking at closely: One can read off that the only possible 
jumps are (x',y') — > (x',y' U y), creation of a boson, and (x',y') — > (x',y' \ y'), an- 
nihilation of a boson. In particular, while one particle is created or annihilated, the 
other particles do not move. The process that we considered in [TT] consists of pieces of 
Bohmian trajectories interrupted by jumps with rates the process is thus an ex- 

ample of the jump rate formula (0), and an example of combining jumps and Bohmian 
motion by means of process additivity. 

The example shows how, for other QFTs, the jump rates (0) can be applied to 
relevant interaction Hamiltonians: If Hj is, in the position representation, a polynomial 
in the creation and annihilation operators, then it possesses a kernel on the relevant 
configuration space. A cut-off (implemented here by smearing out the creation and 
annihilation operators) needs to be introduced to make Hj a well-defined operator on 



If, in some QFT, the particle number operator is not conserved, jumps between the 
sectors of configuration space are inevitable for an equivariant process. And, indeed, 
when Hi does not commute with the particle number operator (as is usually the case), 
jumps can occur that change the number of particles. Often, Hi contains only off- 
diagonal terms with respect to the particle number; then every jump will change the 
particle number. This is precisely what happens in the model of [TTj . 

4 Existence Results 

The configuration space Q is assumed in this paper to be a measurable space, equipped 
with a a-algebra A. Every set we consider is assumed to belong to the appropriate 
cr-algebra: A on Q or the product a-algebra A £g> A on Q x Q. If F is a quadratic form, 
we will usually use the notation (®\F\ty rather than F($, If P(B)ty and P(C)\P 
lie in the form domain of H, we write (V\P(B)HP(C)\V) for (P(B)V\H\P(C)V). 

4.1 Condition for Finite Rates 

For the argument of Section 12.51 to work, it is necessary that (a) the bracket in the 
numerator of (0) exist as a finite signed measure on Q x Q, and (b) the Radon-Nikodym 
derivative of the numerator with respect to the denominator also be well defined. It turns 
out that, given (a), (b) is straightforward. However, contrary to what a superficial 
inspection might suggest, (a) is problematical even when H is bounded. To see this, 
consider the case Jif = L 2 (R) with the natural PVM (corresponding to position) on 
Q = M, and with H the sum of the Fourier transform on Jif and its adjoint, given by 
the kernel 



L 2 . 
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Then, for \l/ real, the bracket in {1} would have to be understood as proportional to 

^(q) cos(gg') *(</) , 

and \l/ G Jf? could be so chosen that this does not define a signed measure onRxl 
because both its positive and negative part have infinite total weight. In fact, \I/ can 
be so chosen that the resulting a( ■ \q') is an infinite measure, o~(Q\q') = oo, for all q', 
and thus does not define a jump process. Note, however, that for \1/ G R £ 2 (IR), 

a( ■ \q') is finite for this H . 

The following theorem provides a condition under which the argument sketched in 
Section |2~51 for the equivariance of the jump rates <r, steps 1-5, can be made rigorous. 

Theorem 1 Let ffl be a Hilbert space, \I/ G M' with = 1, H a self-adjoint operator 
on ffl , Q a standard Borel space? and P a POVM on Q acting on ffl . Suppose that 
for all B C Q, P[B)^> lies in the form domain of H, and there exists a complex measure 
\i on Q x Q such that for all B,C C Q, 

fj,(B x C) = (V\P(B)HP(C)\V) . (64) 

Then the jump rates (JTJ) are well-defined and finite for F-almost every q' , and they are 
equivariant if, in addition, ^ G domain(if). 

Proof. We first show that under the hypotheses of the theorem, the jump rates (P) are 
well-defined and finite. Then we show that they are equivariant. 

To begin with, the measure fi whose existence was assumed in the theorem is con- 
jugate symmetric under the transposition mapping (q,q') tr = (q',q) on Q x Q, i.e., 
/i(A tr ) = fi(A)*. To see this, note that a complex measure on Q x Q is uniquely 
determined by its values on product sets. /i( • tr ) and must thus be the same 

measure since, by the self-adjointness of H, for A = B x C, n(A tT ) = jj,(C x B) — 
(^\P(C)HP(B)\^) = (V\P(B)HP(C)\V)* = fi(A)*. 

We define a signed measure J on Q x Q by J = | Im \i. Let J + be the positive part of 
J (defined by its Hahn- Jordan decomposition, J = J + — see e.g. j21, p. 120]). Since 
fi is a complex measure (and thus assumes only finite values), J has finite positive and 
negative parts. Since /i is conjugate symmetric, J is anti-symmetric. 

We now show that for every B C Q, the measure J(£> x • ) on Q is absolutely 
continuous with respect to P( • ), the "|\1>| 2 " measure defined in (|2*Tj). If C is a P-null set, 
that is (^\P(C)\^) = 0, then P(C)|*) = 0: if P(C) is a projection, this is immediate, 
and if P(C) is just any positive operator, it follows from the spectral theorem — any 
component of ^ orthogonal to the eigenspace of P(C) with eigenvalue zero would lie in 
the positive spectral subspace of P(C) and give a positive contribution to (ty\P(C)\ty). 

7 A standard Borel space is a measurable space isomorphic to a complete separable metric space with 
its Borel cr-algebra. Basically all spaces that arise in practice are in fact standard Borel spaces, and 
so are in particular all spaces that we have in mind for Q (which are countable unions of (separable) 
Riemannian manifolds). Thus, the condition of being a standard Borel space is not much of a restriction. 
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From P(C)V = it follows that (V\P(B)HP(C)\V) = 0, so that J(B x C) = 0, which 
is what we wanted to show. 

Next we show that for every B C Q, the measure J + (Bx ■ ) is absolutely continuous 
with respect to P( ■ ). Suppose again that P(C) = 0. We have that 

S + (B xC)< S + (B xC) + r(B x C) = 
= \I\(B x C) = sup ^ \J(Bi x Cj-)| 

where the sup is taken over all finite partitions {J { Bi = B of B and |J ■ Cj — C of C. 
Now each 1(5; x Cj) = because x • ) « P( • ) and P(Q) < P(C) = 0. Thus 
J+(B x C) = 0. 

It follows from the Radon-Nikodym theorem that for every 5, J + (B x ■ ) possesses 
a density with respect to P(-). The density is unique up to changes on P-null sets, 
and one version of this density is what we will take as a(B\q'). We have to make sure, 
though, that a is a measure in its dependence on B, and from the Radon-Nikodym 
theorem alone we do not obtain additivity in B. For this reason, we utilize a standard 
theorem [23 p. 147] on the existence of regular conditional probabilities, asserting that 
if Q (and thus also Q x Q) is a standard Borel space, then every probability measure v 
on Q x Q possesses regular conditional probabilities, i.e., a function p( ■ \q') on Q with 
values in the probability measures on Q x Q such that for almost every q', p(- \q') is 
concentrated on the set Q x {q'} C Q x Q, and for every A C Q x Q, p(A\q') is a 
measurable function of q' with 

J p{A\q')v{Qxdq') =v{A). (65) 

q'&Q 

We set v{ ■ ) = J + ( • )/J + (Q x Q) and define a as the corresponding regular conditional 
probability times a factor that takes into account that (0) involves the density of J + 
relative to P (rather than to v{Q x ■) or J + (Q x • )): 

a(B\q>) := p(B x Q\q>) '\ q') . (66) 

The last factor exists because we have shown above that J + (Qx ■)<?(•). a( ■ \q') is 
a (finite) measure because p( • |g') is. For fixed B, a(B\q') as a function of g' is a version 
of the Radon-Nikodym derivative d§ + (B x • )/d¥( ■ ) because 

| ( 7(S| g ')PW) 1> / P(B x Q\q')^^l(q')F(dq') = 

q'GC q'&C 

= TOx S ) I p{ B,C W )g^^nBxC). 
q'eQ 
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According to the theorem on regular conditional probabilities that we used, a is defined 
uniquely up to changes on a P-null set of q's. 

Now we check the equivariance of the jump rates a: for any B C Q, 

3? a F(B)^ J a(B\q')¥(dq') - J a(Q\q)F(dq) = S + (B x Q)-J+(Qx B) , 

q'eQ q£B 

using that o is a version of the Radon-Nikodym derivative of J + relative to P. Since JJ 
is anti-symmetric with respect to the permutation mapping (q, q') i— > (q', q) on Q x Q, 
we have that J + (C xB)=J _ (BxC), and therefore 

J^P(S) = J + (5 x Q) - x Q) = S(B x Q) = 

= |lm//(S xQ)i |lm(tf|P(£)if|tf). 

It remains to be shown that F t (B) = (e- iHt/h ^\P(B) e ~ iHt/h ^) is differentiable with 
respect to time at t = and has derivative 



dPt(B) 



dt 



t=o 



Zlm(9\P{B)H\W). (67) 



If \I/ lies in the domain of H, ^ = e lHt ' h \ff is differentiable with respect to £ at £ = 
j2Hl p. 265] and has derivative \i/ = — j^H^. Hence 

~((%\P(B)%) - (V \P(B)%)^ 
=(%\P(B)\(* t - *„)/*) + - *o)/* |P(5)* ) 
converges, as £ — > 0, to 

+ (*|P(S)*> = + i(fl"*|P(S)*> = §Im(#|P(.B)#*}. 

It now follows that J2? CT P = dW/dt, which completes the proof. □ 

We remark that if, as supposed in Theorem [TJ the measure \i exists, it is also unique. 
This follows from the fact, which we have already mentioned, that a (complex) measure 
on Q x Q is uniquely determined by its values on the product sets B x C. 

Another remark concerns how the (existence) assumption of Theorem ^ can be vi- 
olated. Since the example Hamiltonian of Section 13.51 leads to infinite jump rates, it 
also provides an example for which the assumption of Theorem^ is violated, in fact for 
every nonzero \& G J^. To see this directly, note that, while P(B)^> lies indeed in the 
form domain of H (which is Jif since H is bounded), 



(<tf\P(B)HP(C)\ty) = i J dq P-J dq' 



B C 
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where P- J denotes a principal value integral. For B H C = 0, P- J can be replaced by 
a Lebesgue integral. This, together with (|Slj). would leave for /i only one possibility 
(up to addition of a complex measure concentrated on the diagonal {{q,q) '■ q & Q}), 
namely 

But this is not a complex measure for any \l/ since i^*(q) ^(q') /{q — q') is not absolutely 
integrable. This example also nicely illustrates the difference between a complex bi- 
measure v(B, C), i.e., a complex measure in each variable, and a complex measure fi( ■ ) 
on Q x Q: (^/\P(B)HP(C)\^/) is here a complex bi-measure and thus defines a finite- 
valued additive set function on the family of finite unions of product sets BxC C Qx Q, 
which, however, cannot be suitably extended to all sets AC Qx Q. The essential reason 
is that the positive and the negative singularity in l/(q — q') cancel (thanks to the use 
of principal value integrals) for every product set but do not for some nonproduct sets 
such as {(q, q') : q > q'}. In contrast, a (finite) non-negative bi-measure can always be 
extended to a (finite) measure on the product space; see Section EOl 

A related remark on the need for the existence assumption of Theorem One might 
well have imagined that the complex measure fi on Q x Q, extending (j64j) from product 
sets, can always be constructed, at least when H is bounded, as the quantum expected 
value of the bounded-operator- valued measure (BOVM) P x H P on Q x Q, the U H- 
twisted product measure" P(dq)HP(dq') of the POVM P with itself — or, equivalently, 
the product of the POVM P(dq) and the BOVM HP(dq'). Indeed, the nonexistence 
of n for the Hamiltonian in the principal-value example that we have just discussed, as 
well as for the Hamiltonian in the Fourier-transform example at the beginning of this 
section, implies that P x H P does not exist as a BOVM in these cases; if it did, so would 
/i, for all The Fourier-transform example can also easily be adapted to show that the 
product Pi x P 2 of two POVMs need not exist as a BOVM, and in fact does not exist 
when P\ and P2 are the most familiar PVMs for quantum mechanics, corresponding 
respectively to position and momentum. There is, however, an important special case 
for which the product Pi x P 2 of two POVMs does exist, in fact as a POVM, namely 
when Pi and P 2 mutually commute, i.e., when [Pi(B), P 2 (C)} = for all B and C. This 
will be discussed in Section EOI 

4.2 Integral Operators 

In this section we make precise the statement that Hamiltonians with (sufficiently reg- 
ular) kernels lead to finite jump rates. In particular, we specify a set of wave functions, 
depending on H, that lead to finite jump rates. 

4.2.1 Hilbert— Schmidt Operators 

We begin with the simple case in which \l/ is a complex-valued wave function on Q, so 
that the natural configuration POVM P( ■ ) is a "nondegenerate" PVM. What first comes 
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to mind as a class of Hamiltonians possessing a kernel is the class of Hilbert-Schmidt 
operators; for these, the kernels are in fact square-integrable functions on Q x Q. 



Corollary 1 Let Q be a standard Borel space, Jff = L 2 (Q,C,dq) with respect to a o- 
finite nonnegative measure on Q that we simply denote dq, let \l/ G with = 1, let 
H be a self-adjoint operator on Jif , and let P be the natural PVM on Q (multiplication 
by indicator functions) acting on L 2 (Q,C,dq). Suppose that H is a Hilbert-Schmidt 
operator. Then, by virtue of Theorem QJ the jump rates given by are well-defined 
and finite W-almost everywhere, and equivariant. In fact, the jump rates are given by 
(fTUj) with (q\H\q') the kernel function of H . 

Proof. Since if is a Hilbert-Schmidt operator, it possesses an integral kernel K(q, q') 
that is a square-integrable function [2H1 P- 210], i.e., there is a function K G L 2 (Q x 
Q, C, dqdq') such that for all $ G J? , 

HH>(q) = j K(q,q')$(q') dq' . 
Q 

Thus, for all $, $' G Jf, 

($\H\&) = J dq J dq'<f>*{q)K{q,q') = 
Q Q 

(by Fubini's theorem, because the integrand is absolutely integrable) 

dqdq'$*(q)K(q,q') $'(</)• 



QxQ 

It follows that 



(*\P(B)HP(C)\V) = J dqdq'l B (q)y*(q)K(q,q')lc(q')V(q')= (68a) 

dqdq'y*(q)K(q,q')V(q'). (68b) 



BxC 



Note that since H is bounded, its form domain is Jf? and thus contains all P(B)ty. For 
ACQxQ, define 

fi(A) = J **(q)K(q,q')*(q') dqdq'. 

A 

Since 

\y(q)\\K(q,q')\\y(q')\dqdq' 



QxQ 
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n{A) is always finite, and thus a complex measure. (|68J) entails that is satisfied, so 
that Theorem n applies. □ 



We have already remarked that every Hamiltonian H can be approximated by 
Hilbert-Schmidt operators H n . In this context, it is interesting to note that if H is 
itself a Hilbert-Schmidt operator, and if the H n converge to H in the Hilbert-Schmidt 
norm, then the rates a^' Hn converge to a 9,H in the sense that 

J \a*' H -{dq\q') - a*' H (dq\q')\ \V(q')\ 2 dq' ™ 0. 
QxQ 

4.2.2 Complex- Valued Wave Functions 

In addition to the case of Hilbert-Schmidt operators, Theorem applies in many other 
cases, in which the kernel K(q, q') is not square-integrable, nor even a function but 
instead a measure K(dq x dq 1 ). More precisely, K{dq x dq') should be a a -finite complex 
measure, i.e., a product of a complex- valued measurable function Q x Q —>■ C and a 
a-finite nonnegative measure on Q x Q. (Note that this terminology involves a slight 
abuse of language since a cr-finite complex measure need not be a complex measure.) 
The complex measure /i assumed to exist in Theorem is then 

fi(dq x dq') = K{dq x dq') #(c/) . (69) 

This equation suggests that the minimal amount of regularity that we need to assume on 
the kernel of H is that it be a cr-finite complex measure. Otherwise, there would be no 
hope that (J69|) could be a complex measure for a generic wave function that vanishes 
at most on a set of measure 0. The exact conditions that we need for applying Theorem 
Q]to a Hamiltonian H with kernel K(dq x dq') are listed in the following statement: 

Corollary 2 Let Q be a standard Borel space, Jtf? = L 2 (Q,C,dq) with respect to a cr- 
finite nonnegative measure on Q that we simply denote dq, let \1/ e Jrf? with = 1, 
let H be a self-adjoint operator on Jtff, and let P be the natural PVM on Q acting 
on L 2 (Q,C,dq). Suppose that H has a kernel Kidq x dq') for i.e., suppose that 
K{dq x dq') is a a-finite complex measure on Q x Q, and that some everywhere-defined 
version : Q — > C of the almost-everywhere-defined function ^ G L 2 (Q,C,dq) satisfies 

J \^{q)\ \K(dq x dq')\ \^{q')\ < oo (70a) 
QxQ 

P(B)V G form domain(i/) V5CQ (70b) 

(V\P(B)HP(C)\V) = J ^*(q)K(dqx dq')^(q') VB,C C Q. (70c) 

BxC 

Then, by virtue of Theorem Q, the jump rates given by (Q) are well-defined and finite 
F-almost everywhere, and they are equivariant if \1/ G domain(if). 
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Proof. Set 

H(A) = J K(dqx dq')^(q') . (71) 



.4 



The integral exists because of (jT0a|) and defines a complex measure /x, which satisfies 
§U) because of (f7TIcJ) . □ 

We remark that the choice of the everywhere-defined version ^ : Q — > C of the 
almost-everywhere-defined function \l/ G L 2 (Q,C,<ig) does not affect the jump rates, 
since the measure \i is uniquely determined by its values on product sets, which are 
given in ()64|) in terms of the almost-everywhere-defined function \& G M~ '. 

The reader may be surprised that our notion of H having a kernel K seems to 
depend on whereas one may expect that H either has a kernel or does not, inde- 
pendent of ^ . The reason for our putting it this way is that domain questions are very 
delicate for such general kernels, and it is a tricky question for which \l/'s the expres- 
sion (ty\P(B)HP(C)\^) is actually given by the integral ()70c|) . A discussion of domain 
questions would only obscure what is actually relevant for having a situation in which 
Theorem Q applies, which is (J7U|) . Note, though, that if H has kernel K(dq x dq') for 

then it has kernel K also for every ^' from the subspace spanned by P(B)^ for all 
B C Q. 

The conditions (fTUj) become very transparent in the following case: Suppose if is a 
self-adjoint extension of the integral operator K arising from a kernel K(q, dq') that is 
a cr-finite complex measure on Q for every q £ Q and is such that for every B C Q, 
K(q, B) is a measurable function of q. K is defined by 

K$(q)= J K(q,dq')$(q') (72) 
q'eQ 

on the domain Qs containing the $'s satisfying 

J \K(q,dq')\ \&{q')\ < oo for almost every q (73) 

q'&Q 

and 

J K(q, dq') $(</) is an L 2 function of q. (74) 
Q 

That H is an extension of K means that the domain of H contains and HQ = K<& 
for all $ G ^. Then, for a \I/ G ^ satisfying 

y %dg')*(9')^ 2 (2,C,(l?) VBCS (75) 



and 



1 I^Hif^dgOl W)|d9<«3, (76) 



QxQ 
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conditions (J7U|) are satisfied with K {dq x dq') = K(q, dq') dq, and thus Corollary [21 
applies. The jump rates ((TJ) can still be written as in (fTU|) . understood as a measure in 

q- 

Corollary [21 defines a set of good \l/'s, for which the jump rates are finite, for the 
examples of Sections 13.11 13.41 and for ()38a|) . 

4.2.3 Vector- Valued Wave Functions 

We now consider wave functions with spin, i.e., with values in C k . In this case, let 
^*(q) denote, as before, the adjoint spinor, and ^(g) the inner product in C fc . 

CorollaryElremains true if we replace C by C fc everywhere and understand K(dqxdq') as 
matrix- valued, i.e., as the product of a matrix- valued function and a cr-finite nonnegative 
measure. The proof goes through without changes. 

Let us now be a bit more general and allow the value space of the wave function to 
vary with q; we reformulate Corollary [21 for wave functions that are cross-sections of a 
vector bundle E over Q. The kernel is then matrix valued in the sense that (q\H\q') is 
a linear mapping E q i — > E q . 

Corollary 3 Let Q = [J n be an (at most) countable union of (separable) Rieman- 
nian manifolds, and E = [j n E^ the union of vector bundles E^ 1 over Qy 1 ', where 
the fiber spaces E q are endowed with Hermitian inner products, which we denote by 
^(g). Let M' = L 2 (E,dq) be the space of square-integrable (with respect to the 
Riemannian volume measure that we denote dq) cross-sections of E, let \1/ G with 
|| = 1, let H be a self-adjoint operator on ffl , and let P be the natural PVM on Q 
acting on Jtf. Suppose that K{dq x dq'), the product of a a -finite nonnegative measure 
on Qx Q and a section of the bundle [j q , E q <S> E*, over Qx Q, is a kernel of H for ^ ; 
i.e., suppose that some everywhere-defined version \l/ of the almost-everywhere-defined 
cross-section \l/ e L 2 (E,dq) satisfies ()70a|) - (J70c|) (where the integrand on the right hand 
side of (I70c|) should now be understood as involving the Hermitian inner product of E q , 
and ()70a|) as involving the operator norm of K[dq x dq')). Then, by virtue of Theorem 
[TJ the jump rates given by ([TJ) are well-defined and finite F-almost everywhere, and they 
are equivariant if^£ domain(if). 

The proof of Corollary [21 applies here without changes. ()19j) remains valid if suitably 
interpreted. Corollary [31 defines a set of good \l/'s, for which the jump rates are finite, 
for the examples of Sections 13.21 13.31 13.1U1 13.121 and for in case the sum over % is 
always finite. 

4.2.4 POVMs 

We now proceed to the fully general case of an arbitrary POVM. First, we provide two 
important mathematical tools for dealing with POVMs. 

• Any POVM corresponds to a PVM on a larger Hilbert space, according to the 
following theorem of Naimark J7J p. 142]: If P is a POVM on the standard Borel 
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space Q acting on the Hilbert space Jf, then there is a Hilbert space J^ ex t 5 ^ 
and a PVM P cxt on Q acting on J^ xt such that P( ■ ) = P + P ext ( ■ )I with P + : 
^ext - ► the projection and I : Jrf? — > J4? ex t the inclusion, and <¥? cy ± is the closed 
linear hull of {P ext (B)Jf : B C Q}. The pair Jif ext ,P ext is unique in the sense 
that if J$?J xt , Pg Xt is another such pair then there is a unitary isomorphism between 
J#lxt and J^J x t fixing Jif and carrying P ext to P ext ■ 

We call J^xt and P cx t the Naimark extension of Jif and P. We recall that for 
the Hilbert space of positive energy solutions of the Dirac equation and the corre- 
sponding POVM introduced earlier, the Naimark extension is given by L 2 (IR 3 , C 4 ) 
and its natural PVM; this example indicates that the Naimark extension may be, 
in practice, something natural to consider. 

• In Corollaries 121 and El we were considering, for Jri? and P, L 2 spaces with their 
natural PVMs. But when we are given an arbitrary PVM on a Hilbert space, the 
situation is not genuinely more general, since it can be viewed as the natural PVM 
of an L 2 space. We call this the naturalization of the PVM. It is based on the 
following version of the spectral theorem (which can be obtained from the repre- 
sentation theory of abelian operator algebras, see, e.g., jUJ): If P is a PVM on the 
standard Borel space Q acting on the Hilbert space J% ' , then there is a measurable 
field of Hilbert spaces 8 Jif q over Q, a a- finite nonnegative measure dq on Q, and a 
unitary isomorphism U : Jf — > J Jif q dq to the direct integraP of Jtf q that carries 
P to the natural PVM on Q acting on J® J^ q dq. The naturalization is unique 
in the sense that if {JifJ}, (dq)', U' is another such triple, then there is a measur- 
able function f : Q — > (0, oo) such that (dq)' = f(q) dq and a measurable field of 
unitary isomorphisms U q : Jif q — > Jif q ' such that U'^(q) = f(q)~ 1 ^ 2 U q U'^(q). 

A naturalized PVM is similar to a vector bundle in that with every q G Q there 
is associated a value space Jif q , which however may be infinite-dimensional, and 
\I/ G Jif can be understood as a function on Q such that *&(q) G Jif q . Of course, 
instead of the differentiable structure of a vector bundle the naturalization of a 
PVM leads merely to a measurable structure. 

Thus, the situation with a general POVM is not much different from the situation 
with a vector bundle, as treated in Corollary El 

For Hilbert-Schmidt operators, the kernel is so well-behaved that no further condi- 
tions on \l/ are necessary: 



A measurable field of Hilbert spaces on Q is a family of Hilbert spaces Jif q with scalar products 
((■\-)) q , endowed with a measurable structure that can be defined by specifying a family of cross- 
sections such that for all the functions q t— > are measurable and for every q 
the family &i(q) is total in Jff q |18|. 

9 This is the Hilbert space of square-integrable measurable cross-sections of the field {Jtfq}, i.e., cross- 
sections such that all functions q i— ► (($i(q)\$(q))) q are measurable and f(($>(q)\$>(q))) q dq < oo 

im 
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Corollary 4 Let M 3 be a Hilbert space, \I/ G J$? with \\^f\\ — 1, H a self-adjoint operator 
on , Q a standard Borel space, and let P be a POVM on Q acting on M' . Suppose 
that H is a Hilbert- Schmidt operator. Then, by virtue of Theorem the jump rates 
given by (0) are well-defined and finite F-almost everywhere, and they are equivariant. 

Proof. Let P ext be the Naimark extension PVM of P acting on J^ext ^ ^ with P + 
the projection ^f ext — > Jrf?, and let U : Jtf^xt ~^ J%? q dq be a naturalization of P ext . 
For every q G Q, pick an orthonormal basis J? q = {\q,i)} of Jtf'g, with measurable 
dependence on q. When each set J? q is thought of as equipped with the counting 
measure, then from dq we obtain a measure on = M J^, and J J4? q dq is naturally 
identified with L 2 (^,C). Since if is a Hilbert-Schmidt operator, so is H ext = IHP + , 
which thus possesses a kernel function K G L 2 (,y x t y, C) such that for all $ G ^£Lt 



Q 

Since 



C/tf ext <%,i) = f dq' K{q,i,q',i')U$(q',i'). 



(*\P(B)HP(C)\V) = (^\P cxt (B)H cxt P cxt (C)\^) , 
we have, for the same reasons as in the proof of Corollary^ that 

(*\P(B)HP(C)\*) = f dqdq'Y, £ U**(q,i) K(q,i,q',i')U*(q',i') . 

BxC ^ 

For A C Q x Q, set 



(77) 



Since U^*(q, i) K(q, i, q', i') U^(q', i') is absolutely summable and integrable over q, i, q', 
and i', fi(A) is finite, and thus a complex measure. (J77)) entails that (JMj) is satisfied. 
Thus Theorem ^ applies. □ 

We now provide the most general version of our statement about jump rates for 
Hamiltonians with kernel measures. Let t3$(Mqi, J%? q ) denote the space of bounded linear 
operators Jtqi — > Jif q with the operator norm 

\U\ = sup . 

1 /2 

For the norm of \I/(g) in J^, , we also write |\l/(g)|. 

Corollary 5 Let M 3 be a Hilbert space, \& G tyit/i \\^f\\ = 1, H a self-adjoint operator 
on Jtf, Q a standard Borel space, and P a POVM on Q acting on M 1 '. Let P ext be the 
Naimark extension PVM of P acting on J^ ex t 3 and U : Jrf? cxt — > Jtf? q dq the 
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naturalization of P ext . Suppose that H has a kernel K(dq x dq') for \1/ in the position 
representation defined by P; i.e., suppose that K{dq x dq') is the product of a a -finite 
nonnegative measure on Q x Q and a measurable cross-section of the field ^(J^ q ', Jif q ) 
over Q x Q, that \1/ satisfies (|70b|) . and that some everywhere- defined version ^(g) of 
the almost-everywhere-defined cross-section U^/ G J® M'qdq satisfies (j70aj) and ()70c|) 
(where the integrand on the right hand side of ()70c|) is understood as involving the 
inner product of Jtf? q ). Then, by virtue of Theorem^ the jump rates given by (TTJ) are 
well-defined and finite W-almost everywhere, and they are equivariant if^E domain (if). 

The proof of Corollary El applies here without changes if one understands 
$*(q)K(dq x dq')^(q) as meaning ((^(q)\K(dq x dq') $(q'))) q . Corollary El defines a 
set of good \l/'s, for which the jump rates are finite, for the examples of Sections 13.71 
I3~%1 and 1331 

4.3 Global Existence Question 

The rates a t and velocities v t , together with F t , define the process Q t associated with 
H, P, and which can be constructed along the lines of Section 12.31 However, the 
rigorous existence of this process, like the global existence of solutions for an ordinary 
differential equation, is no trivial matter. In order to establish the global existence of 
the process (see [T5] for an example), a variety of aspects must be controlled, including 
the following: (i) One has to show that for a sufficiently large set of initial state vectors, 
the relevant conditions for finiteness of the jump rates, see Sections 14.11 and 14.21 are 
satisfied at all times, (ii) One has to show that there is probability zero that infinitely 
many jumps accumulate in finite time, (iii) One has to show that there is probability 
zero that the process runs into a configuration where a is ill defined (e.g., where the 
denominator of (Tf^|) vanishes, if that equation is appropriate). 

4.4 Extensions of Bi-Measures 

We have pointed out in the next-to-last paragraph of Section 14.11 that a complex bi- 
measure need not possess an extension to a complex measure on the product space, a 
fact relevant to the conditions for finite rates. In this section we show, see Theorem |21 
below, that nonnegative real bi-measures always possess such an extension. 

A useful corollary of Theorem |21 see Corollary [7| below, asserts that one can form the 
tensor product of any two POVMs. This is a special case of the more general statement, 
see Corollary |U] below, asserting that one can form the product of any two POVMs 
that commute with each other; this statement can be regarded as the generalization 
from PVMs to POVMs of the fact that two commuting observables can be measured 
simultaneously; it is also related to the discussion in the last paragraph of Section 14.11 

Though we could not find the explicit statement of Corollary |U1 in this form in the 
literature, it does follow from a part of a proof given by Halmos 20, p. 72]. Below, 
however, we give a somewhat different proof, using Theorem El instead of the lemma of 
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von Neumann [33, p. 167] that Halmos uses. It is also presumably possible to derive 
Corollary El from Lemma 2.1 or Theorem 2.2 of jjj- 

Theorem 2 Let Qi and Q 2 be standard Borel spaces with a-algebras A\ and A 2 , and 
let u( ■ , • ) be a finite nonnegative bi-measure, i.e., a mapping v : Ai x A2 — > [0, a], 
a > 0, that is a measure in each variable when the other variable is a fixed set. Then 
v can be extended to a measure // on Q% X Q2: there exists a unique finite nonnegative 
measure \i : A\ ® A2 — * [0, a] such that for all B 1 G A\ and B 2 G A 2 , 

f i{B 1 xB 2 ) = v{B 1 ,B 2 ). (78) 

Proof. Suppose first that Qi is finite or countably infinite. Then every set A G A\ ® A 2 
is an (at most) countable union of product sets, 

A= |J { qi }xB qi , 
where every B qi G A 2 - Therefore, the unique way of extending v is by setting 

One easily checks that (J79|) indeed defines a finite measure satisfying (J78j) . noting first 
that the sum is always finite because u {{o.i}i B qi ) < Y v({qi}, Q 2 ) = v(Qi, 62)- The 
same argument can of course be applied if Q 2 is finite or countably infinite. 

Suppose now that neither Q\ nor Q 2 is finite or countable. Every uncountable 
standard Borel space Q is isomorphic, as a measurable space, to the space of binary 
sequences {0, 1} N (equipped with the cr-algebra generated by the family B of sets that 
depend on only finitely many terms of the sequence), i.e., there exists a bijection (p : 
<2 — > {0,1 } N that is measurable in both directions, see j2H p. 138] and (231 P- 358]. 
We may thus assume, without loss of generality, that Q\ = {0, l}! -1 - -2 - -3 - -} anc [ Q 2 = 
{0, 1}! ' 1 - 2 '-!, with Bi defined accordingly. Q\ x Q 2 can then be canonically identified 
with {0, 1} Z . 

From the restriction of v to sets B x G B\ and B 2 G B 2 , one easily obtains a consis- 
tent family of finite-dimensional distributions, and hence, by the Kolmogorov extension 
theorem, e.g. [3 p. 24], a unique measure \i on Q\ x Q 2 obeying (|75|l for all B\ G B\ 
and B 2 G B2. 

It remains to establish (|78j) for all Bi G and B2 G ^4.2- First fix i?2 in S2. Then 
/x( • x B 2 ) and i/( • , B 2 ) are measures on A\ that agree on B\. Hence they agree on A\. 
Thus, fixing B\ in A\, we have that /i(-Bi x • ) and v(Bi, ■ ) are measures on A 2 that 
agree on B 2 , and hence on all of A 2 , completing the proof. □ 

In the following, we will again write B\ C Qi instead of B\ G A\. 

Corollary 6 Let ffi be a Hilbert space, Qi and Q 2 standard Borel spaces, and Pi and 
P 2 POVMs on Qi and Q 2 respectively, acting on . If [P 1 (B 1 ),P 2 (B 2 )] = for all 
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B\ C Qi and B 2 C Q 2; ^en t/iere exzsts a unique POVM P on Q\ x Q 2 acting on J*t? 
such that for all Bi C Q x and P2 ^ Q2; 

P(B 1 xB 2 )=P 1 (B 1 )P 2 (B 2 ). (80) 

Proof. (We largely follow |i2TJl P- 72].) For $ 6 we define a bi-measure by setting 
isy(B 1 ,B 2 ) : = (\t / |Pi(-Bi)P 2 (P2) |^)- is obviously a complex bi-measure, and it takes 
values only in the nonnegative reals because Pi(Pi)P 2 (P2) is a positive operator (since 
the two positive operators P\{B\) and P 2 (B 2 ) can be simultaneously diagonalized) . The 
values of v% are bounded by ||^|| 2 . By Lemma |2l can be extended to a measure 
on Qi x Q 2 - 

We now define complex measures on Qi x Q 2 by "polarization": for every 

A <0 Qi x Q 2 and for every pair of vectors $, \& we write 

:= - + »>i*-^*(-A) - (^)- (81) 

We assert that is, for each fixed set A, a symmetric bilinear functional. This 

assertion is proved by noting that (i) it is true if A = B\ x P 2 , and (ii) the class of 
all sets for which it is true is closed under the formation of complements and countable 
unions. To see (ii), note that fx$^(A c ) = /i$,*(Qi x Q 2 ) — ^^(A) and (Ufcli ^k) — 
lim^oo Ylk=i y(A k ). 

Since ^^(A) = fiy(A) < \\^\\ 2 for every AC Q 1 x Q 2 , the bilinear functional 
/j,$ t \s>(A) is bounded and has, in fact, a norm < 1. Therefore, there is a bounded 
operator P{A) such that fi^^(A) = ($\P(A)\iff). P(A) is positive since n^^(A) > 
for every P( • ) is countably additive in the weak operator topology because • ) 
is countably additive. P( • ) satisfies (|HU|). and thus P(Qi x Q 2 ) = J. □ 

Note that Pi need not be a commuting POVM, i.e., possibly [Pi(5 x ), P 1 (C 1 )] 7^ 0, 
and correspondingly for P 2 . 

An immediate consequence of Corollary El which we use in several places of ^2], is 

Corollary 7 Let and Mq, be Hilbert spaces, Q\ and Q 2 standard Borel spaces, and 
Pi and P 2 POVMs on Qi and Q 2 respectively, acting on J^f and ffl^ respectively. Then 
there exists a unique POVM P on Q% x Q 2 acting on M[®J^ 2 such that for all B\ C Qi 
and B 2 C Q 2 , 

P(B 1 x B 2 ) = P 1 {B 1 ) <g> P 2 (P 2 ). (82) 

5 Minimality 

In this section we explain in what sense the minimal jump rates — or (JT§j) or ()38a|) — 
are minimal. In so doing, we will also explain the significance of the quantity J defined 
in (J25|) . and clarify the meaning of the steps taken in Sections 12.41 and 12.51 to arrive at 
the jump rate formulas. 
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Given a Markov process Q t on Q, we define the net probability current j t at time t 
between sets B and B' by 



j t (B, B 1 ) = lim i- [Prob{g t G 5', Q t+At G £?} - (83) 

At\0 At 



-Prob{g t G B, Q t+M G B'} 



This is the amount of probability that flows, per unit time, from B' to B minus the 
amount from B to B' . For a pure jump process, we have that 

3t(B,B')= J * t (B\q') p t (dq') - J a t {B'\q) p t (dq) , (84) 

so that 

j t {B,B')=j a>p {BxB') (85) 
where j^p is the signed measure, on Q x Q, given by the integrand of 

jV,p(d5 x ^') = cr(^k') pW) - o-WIg) p(d?) • (86) 

For minimal jump rates a, defined by (JTJ or (jTHj) or (|38aj) (and with the probabilities p 
given by (JUJ, P = P), this agrees with as was noted earlier, 

jcr,p = , (87) 

where we have made explicit the fact that J is defined in terms of the quantum entities 
\P, H, and P. Note that both J and the net current j are anti-symmetric, J tr = — J and 
j tr = — j, the latter by construction and the former because H is Hermitian. (Here tr 
indicates the action on measures of the transposition (q, q') \— > (q',q) on Q x Q.) The 
property (|H7jl is stronger than the equivariance of the rates a, Jz? CT Pi = dF t /dt: Since, 
by (USD, 

(J? a p)(dq)= Ja , p (dqx Q), (88) 

and, by (J25J) . 

= J(rfg X Q), (89) 

the equivariance of the jump rates a amounts to the condition that the marginals of 
both sides of (j%7|) agree, 

]*A d( l x Q) = Hdq x Q) . (90) 

In other words, what is special about processes with rates satisfying (|%7j) is that not 
only the single-time distribution but also the current is given by a standard quantum 
theoretical expression in terms of H, and P. That is why we call (|87j) the standard- 
current property — defining standard-current rates and standard-current processes. 

Though the standard-current property is stronger than equivariance, it alone does 
not determine the jump rates, as already remarked in 0EH]- This can perhaps be best 
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appreciated as follows: Note that (jHtj)) expresses j u>p as twice the anti-symmetric part of 
the (nonnegative) measure 



C(dq x dq') = a(dq\q') p(dq') (91) 

on Q x Q whose right marginal C(Q x dq') is absolutely continuous with respect to p. 
Conversely, from any such measure C the jump rates a can be recovered by forming the 
Radon-Nikodym derivative 

< d ^ = CJ mr- (92) 

Thus, given p, specifying a is equivalent to specifying such a measure C. 
In terms of C, the standard- current property becomes (with p = P) 

2AntiC = J. (93) 

Since (recalling that J = J + — J~ is anti-symmetric) 

J = 2AntiJ + , (94) 

an obvious solution to (|9*3^) is 

C = J+ 

corresponding to the minimal jump rates. However, (jHTjl fixes only the anti-symmetric 
part of C. The general solution to (|9~3*|) is of the form 

C = J + + S (95) 

where S(dq x dq') is symmetric, since any two solutions to (|9~3|) have the same anti- 
symmetric part, and 5* > 0, since S = C A C tr , because J + A (JI + ) tr = 0. 
In particular, for any standard-current rates, we have that 

C>1\ or <M fl>?&£*Q.. (96) 

Thus, among all jump rates consistent with the standard-current property, one choice, 
distinguished by equality in (jHEJ), has the least frequent jumps, or the smallest amount 
of stochasticity: the minimal rates (Q). 



6 Remarks 
6.1 Symmetries 

Quantum theories, and in particular QFTs, often have important symmetries. To name a 
few examples: space translations, rotations and inversion, time translations and reversal, 
Galilean or Lorentz boosts, global change of phase \l/ — > e l9 ^, and gauge transformations. 
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This gives rise to the question whether the process Q t of the corresponding Bell- 
type QFT respects these symmetries as well. Except for Lorentz invariance, which is 
difficult in that Lorentz boosts fail to map equal-time configurations into equal-time 
configurations, the answer is yes; a discussion is given in jT2 Sec. 6.1]. An essential 
ingredient of this result is the manifest fact that the minimal jump rates ((TJ inherit the 
symmetries of the Hamiltonian (under which the POVM transforms covariantly). 

6.2 Homogeneity of the Rates 

The minimal jump rates (JIJ define a homogeneous function of degree in i.e., cr A * = 
cr* for every A G C \ {0}. This property is noteworthy since it forms the essential 
mathematical basis for a number of desirable properties of theories using such jump rates 
(such as that of [H]): (i) that (when P is a product PVM) unentangled and decoupled 
subsystems behave independently and follow the same laws as the entire system, (ii) that 
"collapsed-away," i.e., sufficiently distant, parts of the wave function do not influence the 
future behaviour of the configuration Q t , (iii) invariance under a global change of phase 
\jz gw^r^ (j v ) invariance under the replacement \& — > e~ lEt / h ^> for some constant E, 
which corresponds to adding E to the total Hamiltonian, (v) invariance under relabeling 
of the particles (which may cause a replacement \& — ► — $ due to the Pauli principle). 

6.3 H + E 

Adding a constant E to the interaction Hamiltonian will not change the jump rates (jTj) 
provided P is a PVM. This is because ($\P(B)EP(C)\^) = E {^\P(B n C)|tf) has 
vanishing imaginary part. For a POVM, however, this need not be true. 

6.4 Nondegenerate Eigenstates 

As mentioned earlier, after (|19|). it is a consequence of the minimal jump rate formula 
(0), in fact of the very minimality, that at each time t either a(q\q') or u{q'\q) is zero. It 
follows that for a time-reversible Hamiltonian H and POVM P, all jump rates vanish 
if \I/ is a nondegenerate eigenstate of H. This is because, in the simplest cases, (q\H\q') 
is real, and the coefficients (q\^) can also be chosen real, or, more generally and more 
to the point, because in this case the process must coincide with its time reverse, which 
implies that the current from q to q' is as large as the one from q' to q, so that minimality 
requires both to vanish. 

6.5 Left or Right Continuity 

From what we have said so far, there remains an ambiguity as to whether Q t at the jump 
times should be the point of departure or the destination, in other words, whether the 
realization t i— > Qt should be chosen to be left or to be right continuous. Although we 
think there is not much physical content to this question, we should point out that de- 
manding either left or right continuity will destroy time-reversal invariance (cf. Section 
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16. ip . A prescription that preserves time- reversal invariance can, however, be devised 
provided the possible jumps can be divided into two classes, A and B, in such a way 
that the time reverse of a class- A jump necessarily belongs to class B and vice versa. 
Then class- A jumps can be chosen left continuous and class- .B as right continuous. An 
example is provided by the model of [TT]: since at every jump the number of particles 
either increases or decreases, the jumps naturally form two classes ("creation" and "an- 
nihilation"), and the time reverse of a creation is an annihilation. The prescription could 
be that if a particle is created (annihilated) at time t, then Q t already (still) contains the 
additional particle. But the opposite rule would be just as consistent with time-reversal 
symmetry, and we can see no compelling reason to prefer one rule over the other. 

7 Conclusions 

We have investigated the possibility of understanding QFT as a theory about moving 
particles, an idea pioneered, in the realm of nonrelativistic quantum mechanics, by de 
Broglie and Bohm. The models proposed by Bell |3J and ourselves |11J turn out to be 
rather universal; that is, their construction can be transferred to a variety of situations, 
involving different Hamiltonians and configuration spaces, and invoking formulas of a 
canonical character. 

One ingredient of the construction is the use of stochastic jumps whose rates are 
determined by the quantum state vector (and the Hamiltonian) . These rates can be 
specified through an explicit formula that has a status similar to the velocity for- 
mula in Bohmian mechanics. We have provided a version of this jump rate formula 
that is more general than any previous one. Indeed, it seems to be the most general 
version possible: we need assume merely that the configuration space Q is a measurable 
space (the weakest notion of "space" available in mathematics), that the Hamiltonian is 
well-defined, and that Q and the Hilbert space are related through a generalized posi- 
tion observable (a positive-operator-valued measure, or POVM, the most general notion 
available in quantum theory of how a vector in Hilbert space may define a probability 
distribution). We have shown that these jump rates are well-defined and finite if the 
interaction Hamiltonian possesses a sufficiently regular kernel in the position represen- 
tation defined by the POVM. 

We have also indicated that in a Bell-type QFT, the different contributions to the 
Hamiltonian correspond to different contributions to the motion of the configuration Q t . 
The relevant fact is process additivity, i.e., that the generator of the Markov process Q t is 
additive in the Hamiltonian. The free process usually consists of continuous trajectories, 
Bohmian or similar, an observation already made in ^1] for the model considered there. 
Exploiting process additivity, we obtain that Q t is piecewise deterministic, the pieces 
being Bohm-type trajectories, interrupted by stochastic jumps. Given a Hamiltonian 
and POVM, our prescription determines the Markov process Qt- As an example, we 
have described the process explicitly for a simple QFT. 

The essential point of this paper is that there is a direct and natural way — a canonical 
way — of devising a Bell-type version of any QFT. 
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